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I. Introduction 

This interdiscipline research program was undertaken in an effort 
to investigate the relationship between the mechanical strength of Mo 
based alloys with their electronic structure. That such a relationship 
may exist for these alloys was suggested by the work of Stephens and 
Witzke.^ The program was to consist of a collaborative effort in which 
Prof. Drew would study the electronic properties of these alloys through 
optical studies and Prof. Arsenault would examine the classical solid 
solution strengthening mechanisms based on size and molecular differences 
to determine if these mechanisms could explain the hardness data. It was 
hoped that these efforts would lead to a theory of the mechanical strength 
of this class of alloys. The program was not successful to this extent. 
However, the optical data does suggest a large potential scattering effect 
in the alloys that increases the more the solute valency differs from 
that of the host. This potential scattering implies a d-+s charge transfer 
which should have a significant effect on the bonding of the solute atoms. 
A more detailed understanding of the mechanical strength of these alloys 
will have to await more theoretical analysis of these results. 

In the next section we describe the optical experiments. Then in 
Section III the work on the classical solid solution theory is discussed. 


2 


II. Optical Experiments 

The best understood disordered transiti'jn metal alloys are those 

2 3 

with noble metal hosts and low concentrations of Ni or Pd. ^ In these 
alloys the impurity d levels fall above the host d bands ^ and lie in 
the nearly free electron s-p bands with which they weakly mix to form 
virtually bound states. The optical properties of these alloys are 
particularly simple as the absorption associated with the impurity d 
levels falls at lower frequencies than the interband edge, and is there- 
fore superimposed on a smooth structureless background. Optical studies 

456 57 8 8—12 8 

of CuNi , ’ ’ AuNi, ’ CuPd, AgPd and AuPd have proven very fruit- 
ful. Also, photoeraission studies on these alloys display the impurity 

13 

density of states clearly separated from the host d bands - 

More generally the impurity d levels fall within the host d bands 

2 

in transition metal alloys. In this case the impurity density of 
states is expected to have a larger energy width and a more complex 
spectrum. The theory of the electronic structure of these alloys is 
much more complex so that our understanding of them is much poorer. Un- 
fortunately the experiments are also more difficult to interpret in this 
case since the host d bands also change upon alloying, and it is difficult 
to separate these changes from those due to the added impurity levels 
(assuming such a separation is justified). Consequently there have been 
few measurements and still fewer interpretation of optical or photo- 
emission spectra for these alloys, even though these are just the alloys 
of the most technological interest. 

In this work we present a study of the optical and electronic prop- 
erties of alloys of the 5d transitional elements with Mo, In particular 
we have studied the MoRe, MoOs and MoPt alloy systems for low 5d element 
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concentrations* These alloys present two kinds of difficulties to optical 
studies* The first is that their electronic structure is complex* The 
second difficulty concerns the preparation of high quality samples of 
these refractory alloys suitable for optical studies* We solved this 
second problem by preparing film by the co-evaporation of the two con- 
stituents of the alloy with separate e-gun sources. The first problem 
has prevented a complete understanding of the electronic structure of 
the alloys but this work has provided important new data which should 
provide significant clues to the determination of the electronic structure 
of these alloys. 

A. EXPERIMENTAL 

The samples were prepared by simultaneous vacuum evaporation of the 
two constituents onto polished fused quartz substrates, 1 mm thick by 
25.4 mm diam. Two substrates were mounted in the evaporation chamber in 
the geometry shown in Fig. 1, so that a shield placed between the two 
sources prevented the solute vapor from falling onto the pure substrate 
while permitting the Mo vapor to fall on both* The deposition rates were 
typically 50 A/sec (monitored by a quartz oscillator) , and the total thick** 

Q 

ness was typically 3000 A. The solute was evaporated with a Varian 2-kilo- 
watt electron-beam evaporation source, while Mo was evaporated! vri.th a 
6-KW electron-beam evaporator. The evaporation chamber was evacuated by 
a conventional oil diffusion pump with liquid nitrogen cold trap- During 
film deposition the pressure was about 5 x 10 ^ Torr and during the sutb- 

^■7 

sequent annealing (400 C for 15 Min) it went down to about 5 x 10 Torr. 
Lower annealing temperatures were unsatisfactory because the reflectivity 
results were then found to change for several days, as the samples con- 
tinued to anneal at room temperature. The 400® C annealing temperature 
produced stable results. The solute concentrations were determined from 
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dc resistivity ratios of the samples between room temperature and 4° K 
using the resistivity data in the literature and calibrated corrections 
by chemical analysis of several of our films. Our concentration determi- 
nation was accurate to about 10% relative. 


Reflectivity measurements were made with a single beam differential 

14 

ref lactometer similar to that described by Beaglehole. Both pure and 

alloy samples were mounted on a rotating holder so that light was reflected 

alternately from a pure molybdenum and an alloy sample. The quantity 

a=(R -R_, )/(R +R-, ) was recorded continuously as a function 

pure alloy pure alloy 

of photon wavelength. 

A typical spectrum of a vs. to is shown in Fig. 2. a was measured to 
one part in 10^ in the visible, and the precision fell off to one part in 

3 

10 in the infrared (below 0.6 eV) and in the ultraviolet (above 5 eV) . 

Only pure and alloy samples prepared at the same time were compared with 
one another in order to reduce the reflectivity differences due to vari- 
ations of vacuum evaporation conditions and the effects of ’aging' after 
exposure to air. The reflectivities of our pure Mo films were measured 
separately using a Perkin-Elmer Model 350 spectrophotometer with a specular 

reflectance accessory. We found agreement to within a few percent over the 

■I 5,16 

measured spectral range with the published results on bulk samples. 

The complex conductivity function a = + ia^ of the sample was ob- 

tained through a Kramers- Rronig analysis of the reflectivity data. O 

X~cx 

and a 1 , were determined separately from R and R = R Tj^T • 

alloy ^ pure alloy pure |i-na; 

The results were then presented in terms of the differential conductivity 

In order to perform the Kramers-Kronig analysis it 


7a. - - 0 ^"’'^ 


1 “1 

was necessary to extrapolate the reflectivity data outside the measured 

spectral range* For the low-energy extrapolation we modeled the complex 

Atria 


dielectric function s = 1 + 


to 


as the sum of a Drude term for the 


conduction electron response and a simple Lorentzian oscillator term to 
account in an approximate way for the interband absorption. The oscillator 
frequency was located well above the low frequency cut off, at 0.5 eV. 
Thus for the low-energy extrapolation, we have 




e(w) ~ 1 - 


0)(£lrt-i/T) 


2 2 
U) - U) 


U < 01 


iFji 


Z 


( 1 ) 


o 

For the Crude term we used the plasma frequency (ui =8.4 eV) given for 

P 
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pure Mo, and assumed that it was unchanged for the alloy. The life-time T 
was taken as the zero-frequency value determined from the dc resistivity. 
The three parameters in the Lorentzian oscillator term were then determined 
by a least squares fit to the reflectivity data near the low frequency cut 
off at =0.5 eV. This procedure is found to provide a physically con- 
sistent low-energy extrapolation. However, the results for Aa are not very 
sensitive to this extrapolation, even if much simpler procedures are used, 
except very near the cut-off frequency. 

We extrapolated our reflectivity data for pure Mo to higher energy 

(oi > 6.8 eV) so that the Kramers-Kronig inverted optical conductivity, 

X6 

a^(o)) = ai matched the data of Koelling et, al. around 

2 eV where their results are most reliable. obtained in 

this manner for our pure Mo film is shown in Fig. 3. For alloy samples 
we modified the extrapolation in order to satisfy the differential con- 
ductivity sum rule 


An = 


2m 


TTe 


Aa^(ui*)do)’ 


( 2 ) 


where An is the difference of the density of valence electrons between pure 


Mo and the alloys. The high frequency cut-off to in Eq. (2) was taken to 

s 
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be the frequency for the saturation of the sum rule for valence electrons 

in the pure system. We chose w to be 30 eV, corresponding to the experi- 

s 
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mental result for pure Au. Although the magnitude of O^((o) within the 
data range is in general strongly dependent on the choice of the high 
energy extrapolation function, the results for obtained 

this way were practically insensitive to the extrapolation functions. 

B. DISCUSSION 

The results for representative MoRe, M 0 O .5 and MoPt films are shown 
in Figures 4, 5 and 6. We will be discussing the following spectral features 
in this section. At low frequencies (below 1 eV) all the curves show a 
sharply rising low-frequency divergence. In the 1-5 eV region there is 
sharp structure in all alloys. At high frequencies (h(0 > 4 eV) the Aa 
spectrum has a sharp rise. The observed spectra were nearly linear over 
the 1-10 at. % concentration range. 

We will discuss the interpretation of the Acr spectra in terms of three 
processes: (a) Modifications in the host density of states which can in- 

clude a reduction due to the replacement of host atoms t*ith the impurity 
atoms and shifts and broadening of host bands, (b) Absorption associated 
with the impurity density of states. This includes transitions from impurity 
states to host states and the modification of the intraband absorption due 
to the impurity. (c) Impurity- induced k-non-conserving interband absorption. 
In this process it is expected that optical edges can occur at frequencies 
corresponding to the separation of peaks in the host density of states and 
the Fermi level, and optical absorption peaks can occur at frequencies 
corresponding to the separation of occupied and unoccupied peaks in the host 
density of states. 
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Of these three processes, (a) and (b) have been considered by Beaglehole 
18 

and Hendrickson in analyzing and interpreting their reflectivity data from 

dilute AuFe alloys. From a theoretical standpoint, only process (b) has 

22 


Several authors have considered 
19,20,21 


been usefully developed in the literature. 

the optical response within the Anderson model of dilute alloys. 

They, have only considered the case where the impurity is placed in a simple 
free electron-like band. Since the final states in our experiment consist 
of both S“p like states and d-like states these results cannot give a complete 
description of our data. For the transitions in which the final states are 
s-p like the result for Acr/c in terms of the Impurity density of states, 


P^CE) is 
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^ I 't 


2„2 

V ^ 2 

CJ 


o 
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-to 


dE Pg(E) 


(3) 


where c is the impurity concentration and s is the electron spin, o) 

P 

is the effective plasma frequency associated with the s-p band and is 
8.4 eV for Mo. to^ is related to the d-s matrix element for the impurity 
and it can be estimated from the strength of the interband edge in pure Au 
to be of order 2-3 eV. 

The low-frequency divergence in Aa can be interpreted large] y in terms 
of the increased scattering of the conduction electrons in the alloy. The 
contribution of the Drude term to the conductivity is 


a = 


1 - itOT 


(4) 


2 ^ 

where O = ne x/m . At the frequencies we are dealing with tox >> 1 and so 
o 


we can approximate (for the sake of discussion) 
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Therefore Aa - ^ A(“) • Eq. 5 clearly gives us a low-frequency diver- 

to 

gence in Aa as observed experimentally. However, we should subtract from 

this a term c(CJ„-o„) where stands for o , and is the Drude part of 
P D P pure D 

the conductivity of the pure sample. When this subtraction is performed from 
Ao calculated from Eq. 3 using A(l/x) estimated from the dc resistivity, we 
get reasonable agreement with the observed low frequency behavior of Ac for 
the alloys. From Figure!! 3, 4 and 5 we see that 1 /t for the alloys gets 
progressively larger as the valence difference between the host and impurity 
increases. This behavior is in direct contradiction to the virtual bound 

I 

state model for alloys as can be seen from Eq. 3. Since the peak in Pg(E) 

should move farther away from the Fermi level for impurities Re, Os and Pt 

p^(0) should decrease respectively. Therefore the first term in Eq. 3 would 
s 

decrease across the series in contrast to the observation. The observed 
behavior is more consistent with the idea that the electron relaxation rate 
will increase as the valence electron number difference (i.e. d-electron 
difference) increases because the impurity than provides a larger perturbing 
potential. 

During the last six months of this research program one of the principle 
investigators (H. D. Drew) visited the University of California at Irvine in 
order to collaborate with Prof. D. Mills on some theoretical questions con- 
cerned with the interpretation of the optical experiments. The result of 

i-.his was an examination of the role of k-nonconserving processes in the alloy 

22 

induced by the addition of the impurities. Qualitatively, these processes 
should lead to the following bands of structures; (1) peaks in Ao corresponding 
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to a transition from a peak in the density of states of the host below 

to a peak above (2) edges in Aa corresponding to transitions from peaks 
r 

in the density of states of the host below to the Fermi level or transitions 

fror.i the Fermi level to peaks in the density of states above E * These struc- 

F 

tures were looked for ^irlth the guidance of band calculations which indicate where the 

peaks in the Mo density of states should lie. The structures in the Aa 

data do not correlate well with the expectations of this model* However, the 

theory of the k-nonconserving process is more complex than indicated above 

and a more definitive evaluation must await realistic calculations based on 

the Mo bond structure* 

The structures in Aa do correlate well with the expectations of 

process (a) above. Since Mo atoms are being replaced by the solute atom 

a proportionate decrease in O is expected as one component of the change in o. 

In the simplest model this contribution would go like Aa -ca * Therefore 

pure 

the peaks in the Mo spectrum at 1*7, 2*3 and 3.8 eV should give rise to 
corresponding dips in Aa. The dips in the observed spectra are in approxi- 
mately the correct positions and, the first two at least, are of approximately 
the correct magnitude. The 3.8 eV peak is anomalous in that it shifts 
steadily to higher energy and grows larger with increased valence difference. 

Now we must discuss the contribution due to the valence electrons of 
the impurity atoms. It is generally expected that the impurity d levels in these 
types of alloys are not sharp virtual levels since they fall within the high 
density of d states of the host* Coherent-potential-approximation (CPA) 
calculations on NiCu alloys and AgPd alloys indicate broad densities of 
states function for the Cu and Ag in the Ni and Pd rich phases* Therefore we 
do not expect sharp features in the Aa spectra from the impurity electrons. 
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Also, from the behavior of 1/t discussed earlier it appears that there is 
significant potential scattering in these alloys which increases with 
valence difference. This potential scattering will account for some of 
the impurity electrons because a s-p scattering phase shift implies an 
s-p like density of states around the impurity atom by the Freidel sum 
rule. 

The MoOs spectrum suggests a broad d-like contribution to Aa located 

around 2-3 eV. This and the -ca contribution would account for the 

pure 

small net Aa in the 1.5-4 eV range. However, this model does not generalize 
for the other alloys. Indeed the large negative peak forMoPt near 4 eV 
suggests another type of interpretation. The 3.8 eV peak in pure Mo results 
from transitions between nearly parallel d-like bands over a large region 
of k-space. The large dip could arise from a shift of the lower band 
toward the Fermi level causing the oscillator strength near 4 eV to 
effectively shift to 2 eV. Such a shift would be consistant with a large 
s-p density at the ?t site in this alloy. 

The MoHe alloy results are more difficult to interpret. The negative 
Aa between 1 and 5 eV is too large to be due to just the component. 

On the other hand the shape of Aa is remarkably close to this explanation. 
There does not seem to be any d-like contribution at all below 5 eV and it 
is hard to understand why that contribution would occur at high photon 
energies. 

In conclusion, we have made a series of careful measurements of this 
differential reflectance between alloys of MoRe, MoOs and MoPt referenced 
to pure Ho in the optical spectral range. The differences in the resulting 
Aa spectra are distinctive but as yet yield no definitive conclusions on 
the electronic structure of these alloys. It is suggested that the s-p 
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density at the impurity site increases with increasing valence difference 

from the host. However, further interpretation must awcit more theoretical 

input* Indeed, these data should provide vital clues to the theorists* 

23 

A preliminary report on the Mo-^e data has been published. The work on 

24 

the Os and Ft alloys was reported an APS meeting. 
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III. Computer Studies of Feierls Stress 

In the Initial research proposal it was proposed that an investigation 
of the Peierls stress be undertaken based on an atomistic computer cal- 
culation. These atomistic calculations were to be based on interatomic 
potentials, and in the case of alloys a different interatomic potential 
was to be assigned to the solute atom. Due to the limitation of the 
support (time and money) it was concluded that the only way any progress 
could be made in these calculations was to obtain someone's existing 
program for making such calculations and modifying it to take into account 
the presence of solute atoms. A program was obtained from Dr. P. Gehlen. 
After much effort we were able to run it on our computer, however, in 
order to do so we had to insert a boundary condition. Later it was 
determined that this boundary condition prevented a determination of 
Peierls strass. 

Concurrent with efforts to obtain the Peierls stress from the inter- 
atomic potential, a computer simulation investigation was undertaken of 
the interaction of an edge and screw dislocation and a size and modulus 
misfit, where the Peierls stress is zero (i.e., represent f.c.c. alloys) and 
also a computer simulation investigation of interaction of size modulus 
misfit interactions with a screw dislocation in a lattice containing a 
large Peierls stress (representing b.c.c. alloys). 

Based on the interaction of the screw dislocations with randomly dis- 
tributed solute atoms, the present results demonstrate that the double- 
kink mechanism would predict both solid solution weakening and strenth- 
ening. In Figure 7, the yield stress at CP K required to overcome the 


original page Ib 
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opposing force due to solute-dislocation Interaction is shovn to be 

lower than that of the double-kink nucleation in the absence of the 

solute atoms. As the concentration of the solute atoms increases, the 

yield stress for the double-kink nucleation decreases. In other words, 

if the rate controlling mechanism of dislocation motion is the nucleation 

of double kink, then solution weakening takes place as a result of 

solute additions. Figure 7 also shows that the stress required for 

the sideward motion of a kink is zero in the absence of solute atoms. 

As the concentration of solute atoms increases, the stress required for 

the sideward motion of a kink also increases. In other xrords, if the 

rate controlling mechanism of dislocation motion is the sideward motion 

of a kink, the solution strengthening takes place by the additions of 

solute atoms. For a given solute concentration, the yield stress first 

decreases to a minimum and then begins to increase for large values of 

E where e represents the size and modules misfit. For a given e , 
t t t 

that is, a given alloy system, the yield stress first decreases to a 
minimum then starts to increase with increasing solute concentration. 
Qualitatively, it is in agreement with the general trend of the experi- 
mental observations. 

It is believed that the strength of the solute-dislocation inter- 
action is determined by the modulus which describes the stress field of 
the dislocation. The modulus for screw and edge dislocations in b.c.c. 
metals, k33, has been expressed in terms of the elastic constants by 
Head. Therefore, k33 instead of G(l) and G(2) for the solute and solvent 
atoms was used to calculate the Interaction parameters for modulus inter- 
action. The elastic constants of the Ta-based alloys have been measured 
using the resonance technique by Armstrong and Mordike. This makes the 
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quantitative comparison for Ta-based alloys possible. The misfit-strains, 
e, of the alloys are calculated from Mordihe. The interaction 
parameters, M, for modulus interactions are calculated from the measured 
value of d’ (G,c)/G(l) as follows: 

For elastic spheres of modulus G(2) imbedded in a matrix of 
modulus G(l), the measured value of d’(G,c) should be; 


iG ^ 2[G(2)-G(1)]G(1) 

dc (l-hi)G(l)+(l-u)G(2) 

v(l) - (1) 

" = 3(l-v(l)) 


From equation 1, the interaction parameter M can be written as 
follows: 


M = 


£(M) 


15(l-v(D) 


“ G(l) 


dc 


( 2 ) 


where v(l) is the Poisson's ratio of the solvent crystal. By substituting 
the corresponding k33 for G(l) and G(2) in the above equations, a comparison 
with experimental results can be made. The experimental data of Ta-based 
alloys by Braithwaite, Mordike and Rogausc.- -r 3.52 solute concentration 
at 77 K is also plotted in Figure 8, where x(p) was given the value of the 
yield stress for the pure Ta at 77 K. Quantitatively, it is in fair agree- 
ment with the numerical results. 

Other experimental results for con^arisons are the Mo-based alloys. 

An investigation was conducted by Stephens and Witzke to determine the 
effects of alloy additions on the hardness of Mo-based alloys. Because 
toere are no adequate data of the shear moduli for the solute atoms, in 
the case of Mo alloys, the shear modulus of the solute and the solvent 
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crystals were derived from the Young’s modulus of the polycrystals by 
using the relationship G(i) = (Young's modulus) /2/ (1+v (i) ) . The experi- 
mental results of micro-hardness tests for Mo-based alloys and the numeri- 
cal results are shown in Figure 8 . Qualitatively, it is in agreement with 
the general trend of the numerical results. Recently the hardness behaviour 
of nineteen binary iron alloy systems was investigated by Stephens and 
Vitzke. Their results indicate that: 

(1) The size effect plays a dominant role in controlling the 
hardness behaviour of binary iron alloys. 

(2) Alloy softening does not occur in binary Fe-Zr and Fe-Hf systems 
because of the very large size difference between solute and solvent atoms, 
that is., the large interaction parameter. This is another evidence 
supporting the general trend of the present numerical results. 
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Figure Captions 


Fig* 1 
Fig, 2 


Evaporation geometry for sumultaneous preparation of pure and 
alloy film of predetermined composition* 


a 


R -R, 
-P A 

R +R, 
P A 


for a Mo-Re alloy. 


Fig. 3 vs for pure Mo. 

Fig. 4 vs ^lOi for a MoRe alloy. 

Fig. 5 ^^^2 ^ ^ MoOs alloy. 

Fig. 6 coAe^ vs for a MoPt alloy. 

Fig* 7 The yield stress at 0® K as a function of total interaction 
parameters. The yield stress was nonnalized to the Peierls 
stress of the base metal. The solid line represents the stress 
required for yielding if only double-kink nucleation is 
rate controlling, and the dashed line represents the stress 
required for yield if only lateral motion is rate controlling. 


Fig. 8 Comparison of experimental results of micro-hardness tests for 
Mo-base alloys with the numerical results. 
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Optical Properties of Mo-Re Alloys’) 

By 

S. K. Bahl and H. D. Drew 

Mfasuirments of tht* reflectivity [/^(purc) — i?(alloy)] of 3000 A thick .Mo-Re films for 
plioton energies between 0.54 and (i eV and Re concentrations around 4 at"„ are rejjorted. 
77 le resistivity of tile film. s is also measured and 2,2 txD cni/at**o found for the elmnjrc 
in resistivity of the films upon alloying. The frequency dependence of the optical dielectric 
function, f = f j — le,, of the pure and alloy films in found from a Kranier.s-Kionig analysis 
of the reflectivity data. The results for the differential optical conductivity c4cji(alloy) - 
- folpiirc)! arc discussed in terms of the changes in the electronic structure of .Mo ujion 
alloying M’ith Re. It is concluded that the rigid band model is inapplicable in this case. The 
results of tlic optical data and the resistivity are consistent with the formation of localised 
Re d'levels 0.8 eV below the Fermi level where pure Mo has a low density of electronic 
.‘States. However, this interpretation is not unique and further studies are suggested. 

Xous presentons des mesures de la reflectivite [J?(pur) — i?(aliage)] pour des couches 
de 30(H) A d‘epaisseur, des energies du photon de 0.54 a 0 eV et une concentration dc Rc 
de 4^\,. Xous avons aussi mesure une difference de resistivite entre Telement puret l‘aliage 
de 2.2 ftO cm/at°„. La dependance en frequence de la function dielectrique optique 
f “ Si — if 2 pour Pelement pur et I’aliage est determinee a partir dhiiie annl 3 ’sc de tj'pc 
Kramers-Kronig des resultats sur la reflectivite. La difference des conduct ivites optique.s 
(^[fifaliagc) — fa(pur)] est rcliee aii.x cliangemcnts de la structure electronique dc !Mo en 
aliagc avec Re. On conclude quo le model a bandes rigides ne s’appliqiie pas dans cc cas. 
I^es nicsures pen vent s’expliquer par 1 ’apparition locale de Rc 0.8 eV cii dcssous du niveau 
dc Fermi au.\ endroits ou ^lo pur n une faible densitc d'etnts elect roniqiies. Toutefois, 
cette interpretation n'cst pas unique et une poursuite de ce.s travenu.x scmhle appropricc. 

1- iiifroductiuii 

The (|ueHtion of the electronic structure of trnn.sition metal alloys is a difficult 
one tliat has nonetheless shown recent progress in both experimental measure- 
mentft and theoretical studies [1]. The series of alloys consisting of Mo })ltis the 
5d transitional elements is an interesting one that has received little attention 
to date. Tlie solubilities are all sufficiently iiigh to permit spectroscopic studies 
and therefore a systematic investigation of tlie electronic structure. Anotlier 
interesting aspect of this alloy series is that their mechanieal properties appear 

10 be related to their eleetronie structure [2]. We have therefore undertaken a 
stmh' of tile ojJtical properties of these alloys. As jinrt of tliis objeetive we have 
carrietl out measurements of tlie optica! properties of Mo- He alloy films, Tlie 
results of these studies are reported in this note, 

2, Experimental Details 

.Mo-!^e alloy films were prepared by e-gun evaporation in a conventional 

011 diffusion ])umped bell jar system having a liquid nitrogen cold trap. The 
base pressure of the .system was typically X 10""Torr which rose to 

Work supported by tlie National At*ronautics and Space .'\<lininistration under (trant 
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2 X 10“®Torr during evaporation* Tlie aouree for evaporation was a boiile of 
Mou,i»oJ^eo,os provided by XASA Lewis Researeli Center* The films were deposited 
on oj)tieally polished fused quartz dist-s (1 in* diameter) preheated to 500 
Tlie evaporation rate was A/s and tlie film thiekncss was typically 3000 A. 
Pure molybdenum films were made by evaj)oration of high-purity Mo (MRC 
Marz grade) using deposition conditions similar to that used for alloy films. 
The composition of the alloy films was determined by wet cliemical analysis 
and by X-ray fluorescence. The structure of the films was studied using X-ray 
diffraction, electron diffraction, and electron microscopy. The films were found 
to be polycrystalline with a lower limit on the crystallite size of oOO A. 

The optical measurements were made using a differential technique, i.e*. by 
comparing the reflectivity of pure Mo film with that of the alloy film. The 
quantity a = + 7/^^) was measured in the energy range of 

0.54 to lieV using an optical system similar to that described by Beaglehole 
[3], Ji and 7^^ are the reflectivities of the pure and the alloy films, respectively* 
For the differential reflectivity measurements samples on 1 in. fused quartz 
discs were used. The spectral dependence of reflectivity of pure Mo films was 
measured using the reflecting stage of a Perkin-Eliner 350 spectrophotometer. 
Due to the geometry of tlie reflecting stage, it was not possible to make ineasure- 
inents on the 1 in. disc samples. Instead Mo films on 2 x 1 in.^ fused quartz 
slides were used* Reflectivities were found to be dependent on substrate depo- 
sition temperatures* Samples prepared at 500 substrate temperature were 
found to have reflectivities essentially identical to the best measurement [4] 
reported on bulk Mo* 

Since the spectral dependences of a and 7?p are experinientalh’ measured, 
this allows us to obtain the reflectivity’ of the alloy as a function of photon 
energy. The resistance of the pure and alloy films was measured at room temper- 
ature and 4.2 K using the four probe method. From these results we deduced 
the change in resistivity upon alloying of Ao 2*2 om/at^t^ for our films. 


3, Results 


The spectral de]»endence of a is shown in Fig. 1. Several films were studied 
w'hose nominal composition were the same. However, variation in the resulting 
amplitude of a indicates that the actual composition varied from film to film. 
For the film corresponding to Fig. I tlie composition was measured, by X-ray 
fluorescence, to be 4 at*\,. The .spectra! features in the optical data for different 
films were found to be identical. We estimate the range of compositions studied 
as 3 to 5 at%* 

The optical constants were obtained from reflectivity data using a Krainer:^- 
Kronig analy.sis which requires apj>ropriate low «0.54 eV) and high (|>fie\') 
energy extrapolations of the reflectivity ont.side the experimental range. 

The ftj 0 extrajioiation was determined by assuming that the complex 
dielectric constant function (t) could be expressed as a sum of Driide term and 
a single harmonic oscillator term whose centre ojq \va*s above the low’ frequency 
cut off Wj = 0*54 eV. Thus for the oj 0 extrapolation* w'e have 


e{(o) 1 + + •* ’■ ’ 

—o)^T -b 10) m,-; — (o- H- iro} 


( 1 ) 


The strength and width parameters / and 7^ were taken as adjustable. The Dnide 
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Fig, 1. Spectral dependence of (i?p i?A>/ 
(i?p 4- of Mo -r 4^0 alloy film de- 
posite on 500 °C fused quartz substrate. 
The thickness of the film is approximately 
3000 A 



Fig. 2, Spectral dependence of w AcoOf the 
film of Fig. 1 as obtained by Kramers- 
Kronig analysis. The arrows mark the 
peaks (-1^) and dip (|) in the conductivity 
for pure Mo 


term in (1) can be calculated knowing the values of plasma frequency (cUp) and 
relaxation time (r) for the electrons. For Mo o>p = 8.4 eV. the value given by 
Veal and Paulikas [4] was used, t was obtained by measuring the residual 
resistivity ratio (RRR) of o at 4*2 to o at 300 K. Assuming tlie conductivity is 
d;^niinated by a free-electron-like portion of Fermi surface, we can write for 
the resistivity, o 


^(300 K) 


7ie-\ri TjJ 


C2) 


wliere T[ and Tj arc contributions to relaxation time due to lattice vibrations 
and iiiipiirities and imperfections in the sample. Other symbols liave their 
usual meaning. At low temperature the contribution to relaxation rate due to 
lattice vibration is expected to be insignificant, Tlie RRR for the pure is there- 
fore 


RRR = . 

+Ti 


T| is determined from the value of dc conductivity [5] of pure Mo wliich can be 
expressed as 


tT(lr 


OJ^Tj 

4 jz 




Using fOp — 8*4 eV, a value for Tj//i of 15.2 (eV)*"^ was obtained. Knowing t|, 
the value of ti and the relaxation time for the film can be determined from the 
RRR of the sample. For the pure film corresponding to the data of Fig. 1 ind 2 
Ave found Tpure = (eV)"^. 

It is assumed for the alloy film that the solute adds another scattering 
mechanism, and then the relaxation rate of the alloy can bo found by a similar 
analysis of its RRR. Ta\]oy was found to be 4.4 (eV)'^. 

From values of o)p and r, the spectral dependence of Drude ierm in (I) AA^as 
calculated. Then the parameters / and F Avere chosen, taking (Oq = 1 eV, to 
match F and di?/dfa at o>i = 0.54 eV for the pure Mo and the alloy films. 
This analysis aa'us performed to provide a consistent loAv-frequeiicy extrapola- 
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tion. The results for Afo= ^2 — ^2 were not very sensitive to tlie 

assumptions made in this analysis- 

The values obtained for parameters /and F were- respeelively, HJS.5 (<'\'/7/)- 
and Lb4 eV/// for Mo, 184 (eV/A)- and 0.(57 eV/A for tiie alloy films. 

For the w -► extrapolation it is assumed that the refleetivity varies as 



for co^ Wh. ^jj is the measured refleetivity at the highest measured f^iergy 

oo 

(o^ = 0 eV. The parameter q was varied until the sum rule ( / oj Afo dro = 

u “ 

= ^TrNe^jm, where A’ is total number of valence eleetrons/eni^) was satisfied. 
For a pure Prude metal above the plasma frequency or for the energy region 
well above the interband transition energy, q is exj)ected to he 4. In our case 
the sum rule required q = 2.552 for Mo and q = 2.522 for tiie alloy. TJiis 
indicated that we had not exhausted all the interband transitions at (3 eV. The 
small difference in value of q of the alloy is due to extra concentration of electrons 
due to the addition of Re. 

The energy dependence of (oAe^ obtained by the Kranier-Kronig analysis 
is shown in Fig. 2. The most prominent feature in the spectnini is the sharp 
peak at about 1.1 eV. However, the smaller structures for energies >1,5 cV 
are significant features as they reproduce for different alloy films. 

4. Discussion 

According to Xoelling et al. [6] the absorption edge observed [4] in Mo at 
=K 1.5 eV involves states at They concluded that if this edge is considered 
in terms of direct transitions (requiring A-conservation) this would involve 
states from Ey to states at higher energ 3 % If, on the other hand, indirect transi- 
tions prevail, transitions from peak in the density of states at 1,5 eA' below L\ 
to states near are to be considered. In terms of the rigid hand model alloying 
Mo with He (containing one extra d-eleetron) should produce the equivalent of 
an upward shift of Ey. In the direct transition model, an iijjward shift of Ey 
should produce a positive peak in (o Aeg at 1.5 eV, and for indirect transitions, 
this shift would lead to a negative peak in lo Afo at 1.5 eV. Since no structure 
in CO A ^2 is observed at 1.5 eV, the reflect 'vity data is therefore inconsistent 
with the rigid band model. 

Another plausible idea for understanding of tlie alloy data is in terms of the 
broadening of the 1.5 eV Mo edge [6] upon alloying. Broadening of the edge 
will give a positive peak below 1.5 eV and a negative peak above 1.5 eV in 
£0 A^ 2 ’ This idea seems to be qualitatively consistent with the data. In fact a 
symmetric positive and negative peak of roughly equal magnitude should be 
expected around 1.5 eV. The smaller values of the negative peak above 1.5 e\^ 
could be due to the interference of this structure with otl: interband structure 

at higher energies. From the change in the slope of the aosorption edge of the 
alloy as compared to Mo it is estimated that the edge is broadened =k 0.1 eV. 
The change in A/t (as obtained from resistivity ratio) is 0.1 eV, a value .similar 
to that obtained for the broadening parameter. This suggests tliat the broadening 
may arise from electron lifetime effects upon alloying. 

On the other hand it aj)pears that the peak in (o Af^ railier too low (i.e. 
^ l.l cV) to be eonsidered part of the edge whieli lias a foot at ^5- 1.24 v\ [4, GJ. 
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Another possible inter]*retation of the data is in terms of transitions fnnu 
localized impurity levels on the Re sites. Thus if the Re d-levels were loeati^d 
about O.S eV below Ky the leading edge of the iieak in to eould eome from 
the onset of transition from this level to states al ove A'j . For a narrow' impurity 
level at we expect [7. 8] 


(0 Afjj 


0(ff> — 


57T(oa Xq 
fo ' 10 * 


(5) 


where 0 is the step function, Xq is the number of impurity slates at and ro,, is 
related to the optical matrix element for tlie transitions from these states. In 
addition w*e w*ould expect a reduced contribution from tlie Mo interhand tran- 
sition because of the reduced Mo concentration in the ulioy. Tlie simplest as- 
sumption for tliis contribution is a term m Afo ^ — cmfolpure)^ i.e.* a uniform 
reduction in io Af^.^ where c in the solute concentration is in al‘\j. Evidence for 
this type of contribution is found in P^ig. 2. Tliere are })eaks in (pure) at 
2.;^ and 4.1 eV and dips at 2.0 and 2.8 eV us indicated on the figure. A good 
anticorrelation of these features wdth the fine struct are in foAso is seen. In tliis mo- 
del, we understand the rapid drop in m A^g near 1,5 v\ as a consequence of tlie 
optical edge in Mo. 

Fur .Vq = 10. ft>a = -i eV, equation (5) predicts <ti — 7 eV at (o ^ 1 eV. 
The experimental value of a Af 2 is eV, u>^\ - *1 eV is close to the value deduc- 
ed [9] from the interband edge in Au and Cii d -> s transitions. On the other 
hand X^^ should be less than 10, As there are only tirl eleetrons in Re having 
ten filled states at an impurity level O.H cV below Ey would lead to an excess 
negative charge on the impurity. However, from the point of viewof gencruli/ed 
Anderson model [10) the impurity density of states is given by 


D{v>) 


r7“b l(n>) 
{c> — Aq — A{(o))- 


AHv>y 




w’herc the width function J(w) — -t ) and is the host electronic 

density of states and Vk\i average of the mixing matrix element. is 

the level shift function [10], If ca — E^ — yl(m) has a root near 0.8 eV below’ A’j.-. 
Avhere Z)p(w) lias a minimum,^ should therefore be small, P'rom (U) w'e see that 
D((o) would therefore be very large near this energy. Moreover, as i^(w) would 
be expected to have weight throughout the s--d band width in Mo we might 
expect some weight above Ey which would enable the satisfying of the cliarge 
neutrality of the impurity. 

The value Ao/r (Ay ^ — ypure) gives some infi>rmation about the 

electronic structure of the alloy which we can illustrate by tlie following over- 
simplified analysis. If we assunie that the conductivity is doiniimted by a free- 
electron-like hand cluiracterized by the plasma frequency measiimi l\v Veal 
and Paulikas wc can WTite [11] 

C I- I 

wlicre 2 is the number of condiudion electrons per atom (foivr)p 8.4 vWz ^ 1), 
ky is the effective P’ernii momentum of the conduction electrons and the are 
tlie phase shifts evaluated at the PVnni level. As.siiniing that the resistivity 
dominated by (>j, because of the d-eharaeter of the impurity atom and taking 
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k*y ^ 1,3 X 10® cm as representative (this is typical of the noble metals) we get 

^ =^ 20 sin- 62 iA£ 2 cni/at^^y . ' J 

This 1ms a maximum value of 20 (lO ein when is i-T/2 corresponding in the = t 

Anderson-Friedel model to a virtual bound state located at the Fermi level. 

iSiiiee tlie measured Ao/r of 2.2 jjtQ em/at is nearly an order of magnitude ; ■ 

sjiialler we take this to indicate that the density of d-states of He at Ey is not 

high. This argument is only suggestive because of the assuniptions made in tlie , v 

aimiysis. However it is consistent with the interpretation of the optical data in 
terms of a loealized level 0.8 e\ below Ey, 

5. Concluding Eemarks 

We see that the differential optical pro^jerties of Mo-Re alloys has two kinds 
of interesting features. Structure above 1.5 eV that seems to indicate an overall 

reduction in tlie.Ho density of states and a large peak near 1 eV that has two ■ 

jjossible interpretations. First, it could arise from localized Re d-states located 

in the valley of the host density of states between tlie bonding and anti-bonding ; 

orbital peaks in thedensity of state. The second possibility is tliat this peak comes 

from a broadening of the 1.5 eV optical edge in Mo upon alloying. For this 

interpretation we are left with the question of where the impurity states are 

located. In particular, the electrical resistivity of the alloys indicates that the 

density of impurity d-!evels at the Fermi level is small. Experiments on other 

alloy.s of Mo with od transition elements would be helpful in distinguishing 

between these different interpretations. We have attempted to study MoOs, 

MoPt. and MoHf alloys by the same technique but without success because tlie 
difference in the evaporation rates for the constituents of the alloy were too 
great to produce useful films. We are modifying our evaporator to permit simul- 
taneous evaporation of the constituents to enable further study of this alloy 
series. 
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Abstract. A new method of extrapolation, in the analysis of Kramers-Kronig relations, 
applying to the case of modulation spectroscopy » is presented. The extrapolation of the 
differentia] reflectance, AR/R, is avoided, since this quantity is not theoretically the simplest, 
and it depends on unmodulated optical constants in an essential way. Rather, the extra- 
polation is made directly in the quantity of greatest theoretical simplicity, the real part of the 
(differential) optical conductivity itself. A formula is found which gives the differential 
optical conductivity in terms of the measured differential reflectivity and the extrapolated 
conductivity. The results are discussed and illustrated by the example of the optical conduc- 
tivity of tungsten, modulated by a chemisorbed hydrogen layer. 


The Kramers-Kronig dispersion relations (KKR) connecting real and imaginary parts 
of the Fourier transform of response functions (susceptibilities) are important in all 
branches of physical science. The KKR express the analytic properties in frequency 
space equivalent to the causal properties of the time-dependent response functions. 

One of the uses of the dispersion relations arises because the susceptibility of greatest 
interest may not be directly measured. Rather, the real part of some analytic function of 
the basic susceptibility is measured. In order to invert the functional relation the full 
complex function is needed, i.e., the imaginary part is needed in addition to the measured 
real part. The KKR is used to compute the imaginary part. 

As an example, the fundamental response function (in the long-wavelength or local 
limit) of a solid, may be taken to be the complex conductivity <r((o) — a'{o>) + i(T"(tu). 
The KKR is 



71 


a'U) 6t/U~ — co^). 
Jo 


Greatest interest attaches to ff'(o)), as this quantity is the easiest one to estimate and 
interpret. However, a' is not measured directly. Rather, the simplest experiments are 
capable of measuring only the reflectivity R, given by R = [r|^ = |(n — l)/(n + l)j^ with 
n = [1 -1- 4;i ia(ct>)/eu]’'^. It is well known that, with r = R* ’■ e‘®, a KKR connecting the 
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real and imaginary parts of In r exists, i.e. a relation determining 6(ctj) from an integral 
over In R. 

In practice a difficulty immediately arises, since R is measured over a finite frequency 
range. It is necessary then to extrapolate In outside the measured range, on the basis of 
theoretical and phenomenological considerations- If the measurements have reached 
asymplotia, or alternatively if one is interested only in sharp structure in tr' which is 
little influenced by the extrapolation, this program can be highly successful. On the other 
hand, there are many cases of interest in which there is a significant contribution to a* 
in the unmeasured region, and where has thus not reached its asymptotic form there. 

On physical arguments, one may have a fair idea of how a behaves in the unmeasured 
region, but toextrapolaie R on the basis of such a guess requires a tedious self-consistency 
procedure which is rarely attempted {D F Aspnes 1976, private communication). It 
would be useful therefore, to have a direct way of combining the measured values of K 
in one frequency range w'ith extrapolated or guessed values of cr" in the remaining range, 
in order to give best values for a* for all frequencies. 

We have produced such a formula for the common case in which the measured quantity 
is a /i?u'ar combination of the real and imaginary parts of cr. This case arises, for example, 
in modulation spectroscopy (Cardona 1 969), where, in effect, th.e small difference between 
the reflectivities of two different samples is measured. Then, calling the difference of 
(effective) conductivities Act we have 

AR/R = 2 Re (Aa/na) = ctAtr' + jSAcr". f 1) 

The so-called Seraphin coefficients (Seraphin and Bottka 1 966), a and /?, may be calculated 
from a knowledge of the conductivity of the unmodulated sample, which is assumed to 
be known, 

The problem which we have posed is: given (by measurement) a, j?, and AR/R over a 
given frequency range L, (symmetrically exter,ded to negative frequencies) and given 
(by extrapolation) Aa* over the remaining frequencies R — L, to find Act' on Z-. This may 
be recognised :is a problem of solving a singular integral equation- The solution is found 
by the methods expounded by Mushkelishvili (1953)- The result is 

Aa’ = 2 Re lim I(ca + i?/) = 2 ReZ^(cu) (2) 

where 


X(z)rr AR/Rit) f Aa'it) 


+ P(r) 


Here 

where 


(Im z # 0). 

(3) 

= (Q(r)/P{r))exp TU) 

( 4 ) 

1 r dt + 

(5) 


and Q, R, as well as P are polynomials. The condition that Z(r) be finite everywhere in the 
complex plane, and the sum rule condition that it vanish as Aaip/Anr as z -* x, where 
AfOp is the modulation of the plasma frequency, determines P, Q, and R sufficiently to 
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give a unique solution, However, in general, for a solution to exist at all requires that 
certain additional sura-rule-like conditions be imposed. 

The general case cannot be categorised succinctly, so we restrict consideration to 
the instance where (f)(a)) = -iln [{a + i/J)/{a — i^)] may be taken as smooth (except 
possibly for a jump at the origin), odd in o>, and less in magnitude than 2n for (oeL, 
At least for the case of modulation spectroscopy it may be shown that (p has the above 
properties. In this case it is possible to choose 0 = R — \,P = 0 without loss of generality, 

A number of points concerning the solution (3) may now be made. 

(i) Since the relation is linear, knowledge of AR/R determines a certain ‘experimentaF 
contribution to Act', which is non-zero only in the interval L. The effect of errors in 
AR/R may be readily discussed. The measured AR/R determines a unique contribution 
to the sum rule for Ac\ 

(ii) The contribution to Atr'to)) for weL coming from the integral over R — L 
may be calculated independently of the measured AR/R, a, and ji. Errors arising from 
poor choice of Act' are easily discussed. The total contribution to the sum rule of the 
extrapolated Aa' is also determined independently. (It is not just A<x' da>.) 

(iii) The integrals encountered in equation (3), although improper, are well defined 
under physically realistic assumptions on Act' and AR/R, X'^(a)) for oj e A — L is real 
and positive as is (oc — iP)X^ for co e L, either diverges or vanishes at the end points 
of L, as a power less in magnitude than unity. If X vanishes, the integrals diverge for r 
at the end point, in such a way that X times the integral is fi nite. I f X diverges, an additional 
sum-rule condition, that the sum of the integrals vanishes when r is at the end point must 
be imposed. In either case, Aa' is continuous at the end points, and is equal to the value 
assigned to it by the extrapolation, 

(iv) There are an infinite number of sum-rule-like conditions which must be satisfied 
by the assumed extrapolation if S{z) and all its derivatives are to be continuous at the 
end points ofL. In our case, these sum rule conditions are not strong and it seems unlikely 
that conditions of continuity on high derivatives than the first will be useful. 

(v) If Aa' can be measured at some point or over some small frequency interval, (say 
by ellipsometry techniques), so that a measured Aa' can be used over part of i? — L, 
it is best to choose the interval in a region where |/?| > \o:\. This gives the maximum 
information. A measurement of Aa' at a point where P cc would give, as a rule, little 
but confirmation that AR/R(x za Aa'. 

In order to illustrate the use of the formula we have applied it to analyse the reflectivity 
of tungsten as modulated by a surface coverage of hydrogen. This experiment has been 
briefly discussed elsewhere (Restorff and Drew 1976) and will be the subject of a further 
report (Restorff 1976). The modulation Ag{oj) has been shown in this case (Fiebelmann 
1976) to be given by 

Aa(cj) = r [anfcu, x) - x)] dx (6) 

^ J-co 

= (2icun/c)aj.(cu) 

where the local conductivity of the hydrogen covered (clean) surface at 

depth X. 

I n figure 1 we presen t the d at a AR /R for the ( 1 00) face as well as a', a" for bulk t u ngslen, 
which is used to calculate a and p. Further discussion of the data* together with data for 
other faces will be published elsewhere. Figure 2 shows the resultant a'(a>) based on two 
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Figure 1. Real and imaginary parts of the conductivity (in frequency units, eV/ft|, for'^^ulk 
tungsten (solid line). Fractional reflectance for the (100) face of tungsten as modulated by a 
monolayer coverage of hydrogen. Experimental error is shown. 



Figure 2. Effective surface conductivity {in frequency units, eV A/ii), of hydrogen covered 
tungsten (defined in the text) as computed from the measured reflectance of figure I and two 
difTcrent extrapolations. Error shown arises from experimental error only, 
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extrapolations* In extrapolation (A), an initial fit using five Lorentz oscillators was made 
to early data including photon energies only up to 5-2 eV. We note that even here the 
extrapolation formula (3) was useful, since a preliminary extrapolation taking cr' — 0 
above 5*2 eV suggests where the oscillators should be centred. (The least-squares fit has 
many minimal solutions which it is not easy to distinguish on the basis of goodness of 
fit.) The oscillator fit above 5*2 eV was then used as the extrapolation in equation (3). 
We also point out that an earlier oscillator fit {Anderson et al 1974), admittedly to some- 
what less precise data, replaced the double peak straddling 5 eV by one broad peak. Note 
also that there is a narrow Lorentz oscillator {not shown) centred at 1 5 eV with consider- 
able oscillator strength, which is necessary to satisfy the sum rule cr' dcu — 35eV^ A. 
This isolated oscillator is an unrealistic feature of extrapolation (A), 

In extrapolation (B) of figure 2, which is based on more recent data up to 5*9 eV, we 
have tried to employ what small amount of theory is available at this point. We first 
note that the new data confirm the splitting at 5 eV and indeed suggest that further 
structure just above 6 eV will be found. 

We do not, however, believe it is meaningful to try to predict structure more than 
about an electron volt above the end of the experimental data. We choose an extrapola- 
tion, cutoff above 20 eV, by = A/(o -f a Lorentz oscillator centered near 6eV, 
The Lorentz oscillator represents the peak which we assume exists around 6 or 6-5 eV 
and A/o) represents a high-frequency average oscillator strength which primarily describes 
transitions from the extra electrons on the hydrogen. The only reason for supposing that 
there is a peak near 6 eV is that all extrapolations show that cr' is rising at 5*9 eV, and it is 
difficult to see how this rise could continue to much higher energies. 

The reason for choosing the form A/cd is that this seems to be typical of asymptotic 
forms for transitions between bound states and the continuum. It could be taken as 
constant v/ithout much change in the results. The reason for choosing the cut-off at about 
20 eV is as follows. 

The wavenumber of the initial and final electron is practically unchanged in the optical 
transition. The energy of plane waves in the vacuum with wavenumbers of order 1/^^ 
is roughly 20 eV above the probable position of the hydrogen level, some 5-10 volts below 
the vacuum. (Plane wave states with k l/a^ in the tungsten have a much lower energy 
relative to the vacuum.) We have assumed that the wavefunction of the electrons on the 
hydrogen siil! have Fourier components with wavenumbers characteristic of the Bohr 
radius, The d-electrons on the tungsten are somewhat more compact, and at the 
surface they are likely to be even more compact than in bulk, according to renormalised 
atom theory. Thus we would expect some contribution at high energies from shifts 
induced in tlie tungsten surface electrons by the absorption of the hydrogen. However, the 
sum rule for this vanishes and since it is expected to be at somewhat higher frequency, 
it should have a comparatively small influence on cr' in the measured region. 

It is not our purpose here to offer a profund analysis of the features to be expected 
in t;'. On the contrary, we suggest that such an analysis is badly needed. What we do 
hope to have shown is how such an analysis can be of use in pinning down the features 
of Afx' over a ‘measured' frequency interval. 

We have made numerical studies, which will be published elsewhere, of various 
extrapolation procedures (see also Balzarotti et al 1975). In this study, we chose a model 
tTj.(o)), from which AR/R was calculated. We added noise artificially, and gave it to the 
analyst to treat as experimental data. (The procedure was 'blind' in the sense that he 
produced a predicted without knowing what the model was.) It was found that if the 
model could be moderately well represented by a few Lorentz oscillators really excellent 
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results were obtained by extrapolation based on oscillator fits (except below ^eV). li 
tails of the type discussed above were used, or other non-Lorentzian shapes employed, 
the *base"line’ of the data was not so good, although many correct features remained. 

Below 5 eV the results were totally unreliable^ The reasons for this are that 
is small and thus relatively uncertain there, and the bulk tr' is large there. Since the extra- 
polation gives contributions roughly proportional to the bulk or' this also introduces 
large errors. 

The example discussed above suffers from the fact that little theory exists on the 
surface conductivities. We also have under analysis the reflectivity difference between 
pure Ni or Mo and dilute NiCu and MoRe alloys respectively. Here we expect to be on 
firmer ground theoretically. One feature in this case is that there will be contributions 
from extra core electrons in the extrapolation. This contribution can be put in with 
ease in the present formulation, in contrast to extrapolations of AJ?/R. 
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I 

Since initiation of the program our efforts have been directed ! 

toward measurements of the differential reflectivity of alloys of moly- 
bdenum and elements in the 5d transitional series. The alloy samples 

have been prepared by evaporation from the bulk in a 10 ^ torr vacuum 
to form opaque thin films on quartz substrate. The difference in re- 
flectivity of the alloy and pure molybdenum similarly prepared is con- 
tinuously recorded as a function of wavelength in the spectral range . 

O O 

from 2500 A to 60,000 A. 

We have measured the differential reflectivity of several samples 
of MoRe and MoOs in the, 5% to 10% concentration range, MoHf has not 
yet been successfully prepared and MoPt ingots have just been obtained. 

We have also measured the reflectivity of our pure Mo films over the 
same spectral range. These results are in good agreement with some 
values reported in the literature but differ somewhat from the recent 
results of Veal .m-'d Poulikar (Phys, Rev. 1974 to be published). We are 
presently studying the Kramers-Kronig inversion of our data for the ex- 
traction of the complex frequency dependent dielectric constant which 
the quantity of direct interest. The most interesting result that is 
already evident in the reflectivity data is the differences in the 
spectra of the MoRe and MoOs spectra. The corresponding differences, in 
the imaginary part of the dielectric constant will hopefully permit us 
to learn something about the impurity d levels in the alloy. In partic- 
ular, the results will test whether the common band model or the virtual 
level model is more appropriate for this alloy series. Data on MoHf and 
MoPt is also desirable for this study and we are presently preparing 
these samples. 


of Ib 
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up to now the program has suffered some personnel problems. 

Dr. Bahl joined the group on May 1 and worked one-half time on the 
project through the summer. Unfortunately Dr. Bahl had to return to 
India for a month and is due to return shortly. Another research 
associate is joining the group in mid October and will be taking over 


the research on a full-time basis 
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During the third quarter of the first year of this program we have 
continued our differential reflectivity measurement?! jn the Mo based alloys 
and worked on the Kramers-Kronig analysis of the optical data* Dr. Bahl 
has returned from India and we are now making rapid progress toward the 
completion of the measurement and anaylsis of the optical properties of 
the dilute alloys* 

In an effort to improve the quality of our evaporated samples we have 
tested samples deposited on a quartz substrate heated to 600° C* This 
procedure is expected to produce a better film because of the outgassing 
of trapped gases and because of the higher mobility of the deposited metal 
atom at the elevated substrate temperature. Having seen a noticeable 
improvement in the reflectivity of pure Mo for samples prepared this way 
we are now preparing all our films on heated substrates. We have completed 
the low concentration measurements of MoRe and MoOs and we are now studying 
MoPt films. 

These optical studies are at the critical stage where interesting 

information concerning the electronic structure of the alloys is becoming 

available. In Fig, 1 and 2 we show the optical data and the preliminary 

results for the imaginary part of the dielectric constant, derived from 

A P 

a Kramers-Kronig analysis of the optical data. ” ^2 ^2 plotted 

in Fig. 2 v. s. photon energy in electron volts. These data show the following 

main features: 1. at low frequencies the increased electronic scattering 

in the alloy produces a positive ft'om the Drude response of the S "baud’* 

electrons near the Fermi level; 2. above about 1 eV the Ag^ becomes negative 

and has the general shape of e 2 of pure Mo. This feature can be interpreted 

interpretation 

in terms of a reduction in the number of electrons in the Mo "d bands". This/ 
raises the question of the location of the impurity d electrons. If they 
go in as relatively localized states (in space and in energy) we would expect 
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them to lie near the ?ermi energy and be partially masked by the Drude response. 

We are refining our analysis of the optical data in order to test these 
hypotheses. The MoPt optical data and the resistivity data which we are also 
measuring for all these alloys, will be very helpful in the analysis. We 
expect significant results from these studies in the very near future. 
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Summary of Yearns Research 


/‘jl€ 


Wo have completed the work on alloy films prepared by evaporating bulk 

2 

alloys. The case of MoRc alloys has been carefully examined. The results of optical 
and electrical measurements indicate that the Re d electrons lie in localized states at 

the Re sice and located about 1 eV below the Fermi level. In addition, it appears that t 
d density of states at the Fermi level decreases upon alloying with Rc, These I 

results can be qualitatively understood in terms of the Anderson extra orbital \ 

1 I 

model as it has been applied to problems in the theory of chemisorption. In 
this theory the impurity density of states can be nonzero over the whole range | 

of the host (Mo) d bands and peaks are expected where the host density of states t 

is low. The Rc d states appear to be lying near the minimum of the Mo density | 

of states near the Fermi level. A short paper on this work has been submitted -h 

for publication in the .Journal of the Less Common Metals and a preprint is 

'1 

attached to this report in the Appendix, ;1 

Evaporating films from bulk alloy samples proved to be only marginally suc- 
cessful. Although it enabled us to complete a study of MoRe alloys, we were 
unable to do a careful study of the concentration dependence of MoRe and it 
proved totally unsati sfactory for MoOs and MoPt alloys. For these last two, 

■0 

the evaporation rates of the two constituents differed too greatly to permit j 

making useful films. Also, difficulties were encountered in ascertaining sample i 

compositions reliably .and inexpensively. 

Because of these difficulties we have installed a second e-gun source in f 

our evaporator .so that wc now have a C> kW source in addition to the original 
2 kW source. With this modification along with other improvements in the evap- ;; 

orator wo can now prepare films of predetermined composition by co-evaporation 
of the constituent metals onto a preheated substrate. The evaporator configur- i 

ation is shown in Fig. 1. With the larger c-gun source higher evaporation rates 
are possible which has produced higher quality films. The pure Mo films are 


now found to have a higher residual resistance ratio 10) and a reflectivity 
that is closer to the best results reported for bulk films in high vacuum. (Fig. 2) 
Examination of these films by x-ray and electron diffraction show that they are 

o 

polycrystalline with crystallite size much greater than 500 A. After installing 
the e-gun and rebuilding the evaporator we have calibrated the film thickness moni- 
tor and determined the optimum conditions for evaporating films. With this improve- 
ment we have solved many of the difficult problems associated with preparing high 
quality alloy samples of known composition for the optical studies. 

We are presently studying MoOs alloys as the logical step following our 
work on the MoRe system. Preliminary optical data for a MoOs film is shown in 
Fig. 3. The strong peak at about 3 eV is indicative of a shift of the solute d 
levels toward lower energies for Os as compared with Re (where the peak is around 
1 eV). Preliminary resistivity measurements for the MoOs films give Ap/c = 1.9 pJ2-cm/at 
Tills value is somewhat lower than for MoRe films for which we have measured 
Ap/c = 2.2 v*ii-cm/at.^. This result is also consistent with the above interpretation. 

We are presently repeating these studies for different Os concentrations and per- 
forming further analysis of these data. 

Considerable interest has developed among the theorists in our Solid State 
Physics group on the question of the effects of the electronic structure of these 
alloys on their hardness properties. Part of this interest derives from the close 
relation between the theories of these alloys and the theories of chemisorption for 
gas atoms on W and As we learn more about the electronic structure of these 

alloys, our effort to explain the dependence of the hardness on the number of d 
electrons in the alloy may be assisted by their participation. 

Wo have almost developed a computer program to determine the Pcierls' 
stress from an anisotropic interatomic potential. Hopefully the inserting of 
a solute atom of a different atomic potential will not be too difficult or too 
time consuiiitng. 

I-roni our classical solid solution modeling wo have been alile to determine 

the effect of a site difference (strain) on the yield stress at 0 K (i-’ig. 4). 

2 

The value of t (Peierls' stress) chosen for Mo was 60 kg/mm . The solid line 


j’L'|)resL’nts tiic yichl stress as a function of size differences tor a solute con- 
centration of 1", whereas, the dashed line represents a solute concentration 

of 10%. The arrows with the chemical symbols on top indicate the strain (size 

difference) of the various solute additions. In the case of Os there is a 

4 5 

difference depending upon the reference employed* ’ 

If a \% solid solution is considered then increasing the size difference 
results in a decrease in strength, i.e., solid solution weakening, which means 
that Hf should produce a larger degree of solid solution weakening than Re. 

u 

This is in direct disagreement with the experimental results, but these 
results were expected, however the modulus effect has yet to be considered. If 
a 101) solid solution is considered then the solid solution weakening occurs 
for strains loss than 0.06, and strengthening occurs for strains greater than 
0.06. Again there is no general agreement between the computer simulation 
results and the experimental results at 0 K. 

We have also determined the strengthening as a function of temperature 
for two strains and two compositions (Tig. 5,6). This higher temperature 
region is of greater interest, for the computer simulation results predict 
strengthening for both values of strain and composition, which is in agree- 
ment with the experimental results. Again, additional compositions will have 
to be examined along with the modulus difference effect. 
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Fig. 1 

Evaporator geometry for simultaneous preparation of 
pure and alloy films of predetermined composition. 
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llpf Ipct Ivi t y of our pure Mo films coinji.ired witli Llio best results 
reporteil for bulk Mu in lil)’li v.icmtm, 'I’lie uc.ret'im'uL ifi seen to be 
excellent bi'low 1 eV .'ind bop.li'ni up, to deviate suljsLanl la I ly only 
at high photon energies. These results are better than any 
previously reported results, on films. 
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Fig* 5 • The yield stress as a function of temperature for two different 

I 

concentrations for a strain of O.OS. 
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concentrations for a strain of 0,09. 
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